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SUMMARY 


An investigation has been made into the ability of a method of integral 
relations to calculate inviscid, 0° angle of attack, radiative heating distri- 
butions over blunt, sonic corner bodies for some representative outer planet 
entry conditions. Comparisons have been made with a more detailed numerical 
method, a time asyn?)totic technique, using the same equilibrium chemistry and 
radiation transport subroutines. An effort to produce a second-order approxi- 
mation (two-strip) method of integral relations code to aid in this investiga- 
tion is also described and a modified two-strip routine is presented. Results 
indicate that the one-strip method of integral relations cannot be used to 
obtain accurate estimates of the radiative heating distribution because of its 
inability to resolve thermal gradients near the wall. The two-strip method can 
sometimes be used to improve these estimates; however, the two-strip method has 
only a small range of conditions over which it will yield significant improve- 
ment over the one-strip method. 


INTRODUCTION 

The method of integral relations is a numerical procedure which has been 
used by many researchers to obtain inviscid shock layer solution around blunt 
bodies in supersonic flow. (See refs. 1 to 6.) The single-strip method runs 
rapidly on the computer and gives results which are within engineering accuracy 
for most applications to adiabatic flow situations. Applications of this 
method to nonadiabatic flow situations involving radiative heat transfer have 
been investigated by Suttles (ref. 1) and Edquist (ref. 2). Suttles studied 
the nonadiabatic flow over two spheres in equilibrium air by using a one-strip 
integral method but there were no detailed comparisons available with other 
numerical methods. Radiative heating contributions were fully coupled with the 
other flow equations in that report. Edquist also used a one-strip method to 
compute equilibrium gas flow over spherically capped cones with sonic corners. 
The radiative heating contributions in that study were not directly coupled to 
the governing flow equations and no detailed comparisons were presented with 
other solution techniques. 

The purpose of this report is to make a detailed comparison of the method 
of integral relations (designated MIR) with a more exact numerical method, the 
time asynqptotic technique (designated TAT) (ref. 7). Radiative heating contri- 
butions are fully coupled with the fluid mechanics in this investigation. The 
merits of MIR as a tool for making radiative heating calculation over blunt, 
sonic corner bodies are evaluated. 

There are two major reasons for an interest in this type of investigation. 
Exploratory missions to the outer planets by instrumented probe vehicles are 
tentatively scheduled for the next decade. The entry environment experienced 
by these probes must be understood in sufficient detail to be able to perform 
atmospheric reconstruction analyses and to insure their survival during a severe 



heating-deceleration pulse. The TAT is currently used for making these calcu- 
lations, but it would be advantageous if a faster running MIR could be imple- 
mented to replace or supplement that method. 

Also, recent investigations into shape changes caused by massive ablation 
during entry into the Jupiter atmosphere have indicated the possibility of pro- 
ducing forebody concavities (ref. 8) which cannot be analyzed by TAT as it is 
now formulated. The ability of MIR to calculate flow over bodies with concavi- 
ties has already been demonstrated by South (ref. 3 ) and Rao (ref. 4 ). However, 
before MIR is used for these studies, the consequences of the cruder approxima- 
tions inherent in the method should be understood. To this end, this study is 
restricted to the solution of problems which are characteristic of planetary 
probes entering hydrogen-helium atmospheres. 

It should also be noted that during the course of this investigation, some 
time was spent in the development of a higher order (two-strip) MIR. A summary 
of the results of this development is presented along with some two-strip MIR 
flow field solutions for inviscid, adiabatic, and nonadiabatic flow over blunt, 
spherically capped cones. 


SYMBOLS 


A,B,C,D,E,F 

a 

Ch>Ch6 

^1 > C 2 , . . . , Cg 

C2’.C3’,C8’,Cg* 

D 

G-j iG2» • • • 1^5 

H 

h 

> ^ 2 » • • • > I 5 

K 

L-j , L 2 , . . . , Lg 
M 


nondimensional constants defined in equations (A 3 ) to (A8) 
speed of sound, a*/V^* ' 

mass fraction of hydrogen and helium, respectively 
nondimensional constants defined in equations (12) 
nondimensional constants defined in equations (17) 
damping factor 

nondimensional functions of p, p, u, v, and H defined in 
equations (2) 

total enthalpy, 

static enthalpy, h*/V^*^ 

nondimensional functions of p, p, u, v, H, 6, and r defined 
in equations (2) 

strip index 

nondimensional functions of p, p, u, v, H, 6, Q, o), 65, r, q, 
and s defined in equations (2) 

Mach number 
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M 

00 

N 

P 

Q 

Qr 

R 

«g* 


s,y 

T« 

u 

V * 

00 

V 

x,z 

“I 

y 

6 

T1 

0b 

0c 

< 

p 

0) 

Superscripts: 

n 

« 


free-stream Mach number, V^*/aoo* 
number of strips 
pressure , 

curvature of reference surface, Rjj*/Rb* 

radiative heat flux, 

radius of curvature, R*/Rfi* 

gas constant, m^/(s^-K) 

radius from axis of symmetry, r*/Rj,j* 

body-oriented coordinates (see fig. 1), s*/Rn*, 

temperature , K 

tangential velocity, u*/V^* 

free-stream velocity, m/s 

normal velocity, v*/V^» 

Cartesian coordinates 

ratio of specific heats 

shock displacement distance, 6*/Rjj* 

transformed y coordinate, y/6 

body angle, deg 

cone angle, deg 

metric coefficient, 1 + Q6 ti 

density, P*/P„* 

shock angle, deg 

iteration number 
dimensional quantity 
undamped quantity 


y»/RN» 


3 



t 


converged value of iterative variable on stagnation line 


Subscripts: 

B 

b 

j 


KO 

K1 

N 

R 

s 

stag 


base 

local body conditions 

refers to particular governing differential equation: 
j = 1 ; shock geometry equation 
j = 2; continuity equation 
j = 3; s-momentum equation 
j = 4 ; y-momentum . equation 
j = 5; energy equation 

refers to conditions on lower boundary (body side) of kth strip 
refers to conditions on upper boundary (shock side) of kth strip 
nose 

radiative 

shock 

stagnation point conditions 
free -stream conditions 


Abbreviations: 

MIR method of integral relations 

TAT time asymptotic technique 


ANALYSIS 

The application of the integral method and the derivation of the approxi- 
mating system of nonlinear ordinary differential equations closely follows the 
development presented in reference 6. However, an Nth order approximation is 
considered herein which involves the construction of N parallel computational 
strips between the body and the shock. This construction of parallel strips 
corresponds to scheme I of the method of integral relations (MIR) as described 
by Belotserkovskiy (ref. 5). For N = 1, the equations which follow will reduce 
to the equations in reference 6. 


Governing Equations 

The conservation equations for the steady axisymmetric flow of an inviscid 
radiating gas written in a body-oriented coordinate system are (from ref. 6) 
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( 1 ) 


_ + + L^ = 0 

• 3s 6 ds 3n 6 3n 


where 


y 

and 6 = 

6(3) 

6 

II = - 
r 

Gi = 0 

Li = -(1 + Q6) tan (o) - 9 t)) 

I2 = pu 

G2 = pv 

L2 = 0 

y 

I3 = p + pu2 

G3 = puv 

d 9 bY 1 

Lo = -I - sin 9 ^ Ip + Qrl4 

\ds ds / 

I4 = puv 

G4 = p + pv 2 

L4 = -(Qr + < cos 9 b)p - Qrpu 2 

I5 = puH 

G5 = pvH 

L 5 = j 


The relationship between the shock and the body geometry in a body-oriented 
coordinate system is governed by equation (1) with j = 1 (ref. 6). 


The relation for L5 has been established by making the local tangent 
slab approximation so that from reference 6, 


3 , 

QRs = — (Qr) = 0 
3 s 

QR = <lRy 

The coordinate system is presented in figure 1 . All these quantities have 
been nondimensionalized as follows: 


y 




u 


u* 


V » 


00 





6« 

Rb» 

%• 

Rb = 

%» 

s* 

%* 

Q = Q»Rn* 



p 


<3R* 


QR = 


P 

00 00 



(V »)2 ^ 

' 00 ' ^ 


( 3 ) 


in addition to equation (1) relationships for p = p(p,h) and for qpy 
as a function of thermodynamic properties across the shock layer are required 
and are presented in later sections. 
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Application of MIR 


In figure 1 the computational field in the s,n coordinate system’ is 
shovm divided into N strips of equal width- Each strip is indexed by an 
integer K where K = 1 for the strip bounding the body and K = N for the 
strip adjacent to the shock. The following linear approximations across each 
individual strip are assumed: 

Ij(n) = (Ij)K1 + [dj)K1 - (Ij)K0](Nn - K) 

p(h) = pki + (pki ” PKoH^n ■ K) 



where 


K - 1 K 

^ n ^ - 

N N 


K = 1,2, ...,N 
j = 2, 3,4,5 

and the Ij are defined in equations (2) . The subscript KO refers to condi- 
tions on the bottom boundary of the Kth strip and K1 refers to conditions on 
the top boundary of the Kth strip. At any given s location it can be shown 
that (Ij)K1 = (Ij)(K+1)0 Pki = P(K+1)0 K + 1 = 2,3,...,N. 

Multiplying equation (1) by 6 and integrating the result with respect 
to h across the Kth strip yields 


K/N 3(Ijr) 


"^(K-1 


)/N 9s 


dn 


d6 pK/N 
ds «J(K-1)/N 9n 


9(Iir) pK/N 

dri 4 ' 


J(K-1 


8(GjKr) 


(K-D/N 9n 


dq 


r 


K/N 


+ ) 

^(K-1 


(K-D/N 


Lj dq = 0 


(j = 2, 3, 4, 5) 


(5) 


Equation (5) is simplified by applying Leibniz's rule to the first term, 
by integrating the second term by parts, and by directly integrating the third 
term. The resulting equation is 


dsV(K-1)/N 
pK/N 
'^(K-D/N 


\ d6 


K/N pK/N 

I-jr dq 

/ ds 

(qljr) 

- \ Ijr dq 

(K-D/N n K-D/N 


(Gj<r) 


K/N 

(K-D/N 


Lj dq = 0 


( 6 ) 



The remaining integrals are evaluated by substituting the linear relation- 
ships from equations (4) and (7) 

r = rij + q6 cos 0b 

K = 1 + Qy = 1 + Q6ri 



After further algebraic manipulation, the transformed governing equations 
are written: 

Continuity: 


Ci^(Pu)ko - (Pu)kiJ + C2(pv)ki + C3(pv)kq + C4(pu)kq + C5(pu)ki 

d(pu){{0 d(pu)Kl 

+ Cfi + Cy = 0 

ds ' ds 

Normal momentum: 


Ci|j(puv)Ko - (Puv)ki] + C2(p + Pv2)ri + 

*d(puv)Ko 


+ C5(Puv)ki + C6 


ds 


- Q(p + pu2)|fo 


C3(p + pv2)j{o + C1 |(puv)ko 


■ 

'd(pUv)Ki ' 

+ Cy 

- Q(p + Pu 2 )ki 
ds 


+ C8PK0 + CgPKi = 0 


( 8 ) 


(9) 


Tangential momentum: 

C-|[(p + pu2)fjo - (p + pu2)jj-|] + C2 (puv)ki + C 3 (puv)ro + Cn(pu2)jjo + C5(pu2)jj^ 



"d(p + Pu2)ko 


"d(p + pu2)jfi 

+ C6 

+ Q(puv)ko 

ds 

+ Cy 

+ Q(puv)xi 

ds 


Energy : 


Ci^(PuH)ko - (PuH)kiJ + C2(pvH + qR)Ki + C3(pvH + qR)K0 + Ch(puH)ko + C5(puH)ki 
d(puH)K0 d(puH)Ki 

+ C6 + Cy = 0 (11) 

ds ds 
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where 


Cl = [3N(2K - Drij + 2(3R2 - 3K + 1)6 cos 9b]“ 
C2 = 6n2^1 + Q 6 + 6 | cos 9bj 

C3 = -6n2^1 + Q 6 


K - 1 ) 

+ 6 cos 0K 

N 


C4 = 6 


drb d0b 

3N (3K - 2)6 sin 0b 


ds 


ds 


drb d0b 

3N (3K - 1)6 sin 0b 

ds ds 


C5 = 6 

G6 = 6[3Nrb + (3K - 2)6 cos 0b] 




( 12 ) 


C7 = 6[3Nrb + (3K - 1)6 cos 0b] 
C8 = -6 cos 0b[3N + (3K - 2)6Q] 
Cg = -6 cos 0b[3N + (3K - 1)6Q] 


It has been noted (ref. 6) that numerical instabilities can be elimi- 
nated for N = 1 if the integrated form of the tangential momentum equation 
(eq. (10)) is replaced by the exact form of the tangential momentum equation 
on the body. 


The exact form of the tangential momentum equation on the body can be 
obtained from combining equation (1) (j = 2) with equation (1) (j = 3) while 
noting that q = 0 and v ,= 0 to yield 


^(13)1,0 

ds 


ui,0 


d(l2)i 0 

— = 0 

ds 


( 13 ) 


When K = 1, equation (13) replaces equation (10). 
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Solution Procedure 


Overview .- The ordinary differential equations (eqs. (8) to (10) and (13)) 
are numerically integrated around the body from the axis of symmetry (for 0° 
angle of attack) , adiabatic flow being assumed (H = H^) . The limiting forms of 
the governing equations are obtained on the axis of symmetry and it is found 
that for N strips, there are N more unknowns than equations. The N addi- 
tional relations required on the axis of symmetry are obtained from regularity 
conditions on the sonic line and a sonic corner condition which must be satis- 
fied downstream. Consequently, the solution procedure is similar to a shooting 
method for a two-point boundary value problem, in which N initial conditions 
on the axis of symmetry must be iteratively adjusted until N - 1 regularity 
conditions and a sonic corner condition (Miq = ‘'•0) ai'e satisfied downstream 
on the sonic line. 

For N = 1 , the shock standoff distance ^stag initial condition 

which is iteratively adjusted to force M = 1.0 at the sonic corner. For 
N > 1 , the velocities at the boundaries of adjacent strips are also iteratively 
adjusted so that the governing equations can be integrated across the saddle 
point singularity located on the sonic line. (A description of the saddle point 
singularity eind the behavior of the governing equations near the singularity 
are presented in refs. 3, 5, and 6.) 

A fully converged adiabatic solution is obtained when all downstream regu- 
larity conditions and the sonic corner condition are satisfied by using the 
adjusted iterative variables at 3=0 for initial conditions. The thermody- 
namic conditions around the body are saved and radiative heat fluxes are calcu- 
lated and saved for use in the energy equation. Equations (8) to (11) and (13) 
are now numerically integrated around the body from the axis of symmetry (by 
using the updated values of the radiative heat fluxes) and a new set of initial 
conditions must be obtained. When a converged solution is obtained by using 
the old values of radiative heat fluxes, new values for the radiative heat 
fluxes can be calculated by using the updated thermodynamic distributions. 

Thus, there are two iterative procedures involved in obtaining a converged non- 
adiabatic solution. The first procedure (Procedure I) involves an iterative 
process on the initial conditions to obtain a solution around the body while 
holding the radiative heat flux distribution unchanged (even though thermody- 
namic distributions are changing during the process). The second procedure 
(Procedure II) involves the process whereby the radiative heat fluxes are 
updated everywhere in the field after which program control is passed back 
to Procedure I. As program control alternates between Procedure I and Proce- 
dure II, changes in the radiative heating distribution go to zero and conver- 
gence is achieved. 

Initial conditions .- The transformed system of governing equations includes 
4N ordinary differential equations (eqs. (8) to (11) and (13)), which are inde- 
terminant at s = 0. The limiting forms of these equations are obtained by dif- 
ferentiating them with respect to s. However, all the terms in the tangential 
momentum equations vanish quadratically at s = 0; thus, a single differentia- 
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tion does not remove the indeterminacy. One can either differentiate the tan- 
gential momentum equations again with respect to s and solve for d^p/ds^ 
along the symmetry line on the strip boundaries or ignore the tangential momen- 
tum equations on the symmetry line. The first approach ultimately requires a 
specification of d^p/ds^ on the stagnation point whereas the second approach 
only requires a specification of p on the stagnation point. Because d^p/ds^ 
could not be adequately specified for many different body shapes, the second 
approach (ignoring tangential momentum on the symmetry line) was implemented. 

The limiting form of the governing equations on the symmetry line are: 
Continuity: 


C2’(pv)xi + C3 *(pv)ko + 2C4jp — ) + 2C5[p — ) = 0 

\ ds/Ko \ ds/Ki 


Normal momentum: 


(14) 


C2'(p + PV^)ki + C3'(p + pv2 
- C5QPKI + Ca'PKO + Cg'pKi = 0 


( du\ / du 

)kO + 2C4 pv — + 2C5 pv — 

\ ds/KO \ ds^ 


■) - C4QPKO 

/K1 


(15) 


Energy : 

C2'(pvH + qR)K1 + C3'(pvH + qR)K0 + 2C4(pH — ) + 2C5(pH — ) 

\ ds/KO \ ds/Ki 


= 0 


where 


C2’ = + Q^stag "j(~ - '^stag — 


ds N/ 


ds N 


■ p/ K - l\/drb dSb K - r 

03' = -6N 1^1 + Q6gtag - ^stag ~ 

d0b 

C8’ = «stag ^3N + (3K - 2)6gtagQ] 
d0b 

Cg ' = Sstag + (3K - DSstagQj 


( 16 ) 


(17) 
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The unknovms in equations (14) to (16) are 


Unknown 

Number of unknowns 

Pio» P20 Pno» Pni 

N + 1 

V10» V20. •••» VNO. VN1 

N + 1 

duio/ds, dU2o/ds, ..., dufjo/ds, du(ji/ds 

N + 1 

QRIO 1 QR20» •••> QRNO. <5RN1 

N + 1 

H101 H20 Hno, %i 

N + 1 

Pio> P20> •••» pnoi Pni 

N + 1 

^stag 

1 


The variable is equal to from conservation of energy across a 

shock. (For adiabatic flow H is a constant equal to and equation (16) 

is redundant.) Values of qRio> QR20> •••» <1RN0> QRN1 equal to zero 

for an adiabatic solution and are calculated from thermodynamic data on the 
stagnation line (obtained from Procedure I) by using a tangent slab approxima- 
tion for nonadiabatic solutions. Note also that v-|o = 0 (no mass flux through 
body) and that the stagnation pressure p-]Q can be approximated as (ref. 6) 

1 2 
PIO = Ps + - PsV 

In addition, there are N + 1 local equations of state (see appendix) 
relating pressure to density and enthalpy on the strip boundaries. There are 
also two more equations (Rankine-Hugoniot relations) governing the conservation 
of mass and momentum across a normal shock. Consequently, there are a total of 
4N + 3 equations arid 5N + 3 unknowns. The N additional relations required 
come from regularity conditions on N - 1 strip boundaries between the body 
and the shock which must be enforced on the sonic line and the imposition of 
a Mach number of 1.0 at the aft comer of the body . There are N unknowns 
^^stag» ^20f ^30> •••» which must be guessed at initially and 

iteratively corrected until all N conditions are satisfied simultaneously 
(Procedure I). 

This procedure is most easily implemented for N = 1 where, in the sonic 
corner problem, only one variable is iterated to force a Mach 1.0 condition at 
the aft corner of the body. If the Mach number M-|o is greater than 1.0 
before the corner or if the solution to the governing equations will not con- 
verge at some downstream location (see next section), 6stag increased. 

If the Mach number is less than 1.0 at the aft corner of the body, Sstag 
decreased . 
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Solution of governing equations downstream of symmetry line .- The deriva- 
tives in equations (8) to (11) and (13) are calculated by using three-point 
backward finite-difference formulas. For example, 

= CiPu-|o(s) + C2Pu-|o(s - ^s-]) + C3Puio(s - (^3-| - As2) 

10 

where 

Cl = (1 + 2f)/[f As 2(1 + f)] 

C2 = -(1 + f)/(f As2) 

C3 = y[As2(1 + f)] 
f = Asi/As2 

Asi current increment in s 
As 2 previous increment in s 

The conservation equations, coupled with the local equations of state (appen- 
dix A) , the Rankine-Hugoniot relations (ref. 6) , and the shock geometry rela- 
tion (from eq. (1), j = 1) 

d6 

(1 + Q6) tan (o) - 65) = 0 ( 18 ) 

ds 

now form a system of 5 N + 5 nonlinear, algebraic equations with 5 N + 5 
unknowns. These algebraic equations au’e solved at every integration step by 
using a quadratically convergent Newton-Raphson technique. For the case of 
adiabatic flow, H is a known constant and the energy equations can be 
replaced; thus, a system of 4 N + 5 equations with 4 N + 5 unknowns remains. 
The equations cire solved at a given s location, an increment As is chosen, 
and the procedure is repeated at a new body station until the sonic corner is 
reached or the Newton-Raphson technique fails. (This failure occurs because 
the system of equations becomes ill-conditioned.) At this point the initial 
conditions eire adjusted and the entire solution procedure is repeated until all 
the regularity conditions on the sonic line are met and the Mach number at the 
aft corner of the body is' equal to 1.0. 

Radiative heating calculations .- The coupling of the radiative heat flux in 
the energy equation is accomplished in an iterative manner (Procedure II). An 
adiabatic solution is obtained and the thermodynamic state at each computational 
point is saved. The radiative heat flux along a ray normal to the body is cal- 
culated from these data by using the tangent slab approximation as utilized in 
subroutine RAD/EQUIL (ref. 9 ). Within RAD/EQUIL, line and continuum radiation 
are included. The radiative heat fluxes on the strip boundaries are stored and 
used in the energy equation to compute an approximation to the nonadiabatic 
solution. A new thermodynamic state from the latest nonadiabatic solution is 


dpu 

ds 
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used to compute new radiative heat flux distributions. This successive substi- 
tution process is continued until 


n n-1 

^R,stag “ ^R,stag 


< 


n 

^RjStag 


0.0005 


Values of qp away from the stagnation point can have larger differences 
between iterations but it was found that the convergence criteria specified at 
the stagnation point will result in all variables being converged to within 
three significant figures . A flow chart of the entire solution procedure is 
presented in figure 2. 

In some cases with strong radiative heating, it was found necessary to 
damp the calculated radiative heat flux when substituting into the energy equa- 
tion. Thus, it would be specified that 

Qr” = Dqp'^ + (1 - D)qR'^“^ 

where D = 0.75 would usually be a sufficient damping factor. For all the 
cases considered, the convergence criterion for the nonadiabatic solution was 
usually satisfied after five to seven successive substitutions. 


Two-Strip Solutions 

From a practical point of view, MIR is limited to a one -strip solution 
(N = 1). The complexity of the iteration procedure makes multistrip solutions 
very difficult to obtain. The problems of a multistrip approach to the method 
of integral relations can be summarized as follows by using a two-strip approx- 
imation with equally spaced strips as an example. 

In the two-strip approximation, there are two initial conditions on the 
X stagnation line which must be specified: 5stag» shock standoff distance, 

and V20» the velocity on the strip boundary between the body and the shock. 

A solution is obtained when the governing equations can be integrated along 
each strip to a line normal to the aft corner of the body with a Mach number 
of 1.0 at the corner. For the conditions illustrated in figure 3i the equa- 
tions must be integrated across the sonic line. The position of the sonic 
line in the flow field gives an indication of the location of singular points 
of the governing equations (ref. 5). The iterative variable <Sg^^g and V 20 
must be adjusted to a precision fine enough to simultaneously allow the gov- 
erning equations to be integrated through a singularity in the transonic region 
of the flow field and to force the Mach number condition at the corner. 

In reference 6, it was reported that about four figure accuracy in ^stag 
is necessary to integrate through the transonic region of a sphere using a one- 
strip MIR. It was also reported that there was poor agreement with experimental 
data in the transonic and low supersonic regions because no criteria can be 
established to guide further refinement of 5 stag improve the solution. In 
the flow field depicted in figure 3 , much of the region above the cone flank 
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can be classified as transonic and low supersonic. Therefore, in a two-strip 
procedure, one should expect some difficulty in using MIR since the upper strip 
will pass through a long region of transonic and low supersonic flow. 

Many of the difficulties associated with a two-strip MIR can be overcome 
if the flow field surrounding the body is characterized by a sonic line which 
lies very close to the downstream boundary of the computational domain (fig. 4). 
Such flows can be found on very wide angle bodies although the details of the 
location of the sonic line are also Mach number and y dependent. A two-strip 
solution procedure with equally spaced strips was successfully implemented for 
flows exhibiting this type of behavior in helium (y = 1.67). The criteria for 
obtaining a solution can be summed up as follows: 

(1) M-|o = 1 ± 0*01 on aft corner of body 

(2) Shock angle to is monotonically decreasing 

(3) V2q(s) < 0 

The iterative variable 6g^.gg is held constant^and the ^variable V20 is 

V20 - V20 

iterated according to the criteria described. If < 10“°, where the 

V20 

superscript n indicates the nth iteration, and all of the solution criteria 
have not been met, then 6gtag is adjusted and the iteration on V20 is 
implemented again. The procedure is continued until all the solution criteria 
are met . A sample calculation is presented in figure 5 with comparisons to 
experimental data from reference 10 that use this approach. 

However, there are still some problems with this approach. The execution 
time for this program on the computer is significantly increased because it 
involves a double iteration procedure cind because a smaller integration step 
size (as compared with the one-strip MIR) is required to maintain stability. 

The solution procedure is not closed. There is a small range of values for 
*^stag ^20 which can be used that satisfy the solution criteria. In the 

small number of cases tested, this relationship between Sgtag ^20 

linear (fig. 6). The final problem is that the iteration procedure will only 
work for values of 6gtag ^20 dose to the true solution. A one-strip 

solution must be obtained to determine an approximate value for Sgtag I'® 
used in the two-strip solution. Then a range of values for V20 must be 
determined so that the solution can be integrated far enough downstream to 
insure that the criteria for adjustments on V20 are still valid. 

Typically, values for V20 differing by 5 percent from V20 could be 
iterated to obtain a solution to the problem. (The quantities V20 and 
6stag ®^^® 6he values of V20 and 6stag which satisfy all the convergence 
criteria.) Values of V20 differing by more than +5 percent from V20 
would cause the solution to diverge before the sphere-cone junction and no 
clear reason for this divergence was observed. In many cases, this divergence 
was manifested by the fact that the Newton-Raphson solution of the governing 
system of nonlinear algebraic equations would fail to converge to any solution 
near the stagnation line. (Other two-strip solutions have been reported by 
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Belotserkovskiy (ref. 5), Inouye and Marvin (ref. 11), and Holt and Hoffman 
(ref. 12). The complexity of the iteration procedure and the increased com- 
puter time required for mathematically closed two-strip solutions have caused 
most researchers using the method of integral relations to concentrate on 
obtaining as much information as possible by using the relatively simple one- 
strip approach.) 


For wide-angle bodies where the sonic line lies close to or outside of the 
computational plane (fig. 3), a workable two-strip procedure can be created by 
eliminating the requirement of equally spaced strips. When a very narrow strip, 
Aril = 0.075, is placed close to the body and a second wide strip is used to 
complete the shock layer , it is found that solutions can be obtained by using 
the same integration step size as the one-strip methods. Furthermore, the 
solution became relatively insensitive to the choice of one initial condition 
V20 (the value of V20 the stagnation line which satisfies all conver- 
gence criteria to the sonic corner). Consequently, the solution procedure 
behaves like a one-strip MIH if a reasonable specification can be made for 
V20* assuming a quadratic pressure variation across the stagnation line 
3p 

with — = 0 at ri = 0, one can write 

3n 

P20 ~ Pstag + (Ps “ Pstag^(^Hl)^ (19) 

This value of P20 can be used to calculate V2o from the limiting form of 
the governing equations on the stagnation line. The reformulation of this two- 
strip problem requires another derivation of the governing equations. Now, 
rather than integrating from (K - 1)/N to K/N in equation (5), it is neces- 
sary to substitute the proper values of q at each strip boundary. The only 
- change vdiich results is in the strip constants of_ equations (12). These strip 
constants are now expressed as follows: 


Ci = BPA2 rt, + BB3 6 cos 65 

C2 = (1 + 6Q B1)(rt, + 6 B1 cos 9b) 

C3 = -(1 + 6Q A1)(rb + 5 A1 cos 85) 

/ dTb d8b\ 

Cn = 6 BMA2 - BB1 6 sin 9b 

\ ds ds / 

f drb d9b\ 

Cc = <S BMA2 - BB2 6 sin 9h 

\ ds ds / 

C6 = 6(BMA2 rb + BBT 6 cos 9b) 

C7 = 6(BMA2 rb + ^ 6 cos 9b) 

Cg = -6 cos 9b(BMA2 + BBT 6Q) 

C9 = -6 cos 9b(BMA2 + BK 6Q) 



( 20 ) 
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where 


BBT = (bT + 2rT)BMA2/3 
BB2 = + aT)BMA2/3 

w = (bT 2 + rr bt + rr2)/3 

BPA2 = (bT + aT)/2 
BMA2 = (bT - aT)/2 

rr = (K - datii 

bT = Am + (K - 1)(1 - Am) 

The two-strip solution procedure described requires adjustment on only one 
initial condition and has been used as a tool for assessing the effectiveness 
of the MIR for use in predicting radiative heating over sonic corner bodies. 


RESULTS AND- DISCUSSION 

Flow field calculations have been made over spherically capped cones and 
hyperbolic bodies of revolution with the numerical techniques described herein. 
All the body shapes considered are characterized by having the sonic line of 
the flow field attached at the aft corner of the body. All the solutions are 
obtained for hydrogen -helium gas mixtures because of the current interest in 
exploration of the outer planets by instrumented probes. 

A time asymptotic technique (TAT) as developed by Sutton (ref. 7), for 
inviscid, radiating flow fields has been used as a basis for all the compari- 
sons. The technique, which has been tested against other flow field routines 
and experimental data, has been found to be an accurate, reliable method which 
has been used extensively for making radiative heating calculations over plane- 
tary probes (ref. 13). The TAT was initially developed for treating smooth 
bodies without sonic corners, but a program modification allows a special treat- 
ment of the governing equations at the aft corner of the body to treat a sonic 
condition by following the method of Barnwell (refs. 14 and 15). However, the 
program modification fails when the sonic line of the flow field moves outside 
of the computational domain of a body-oriented coordinate system. This failure 
occurs because the backward difference formulas used to obtain the outflow 
boundary conditions in the TAT are only valid in a supersonic domain. (Barnwell 
uses polar coordinates at the corner to overcome this problem, but Sutton's 
program was never modified in this manner.) Consequently, for given free- 
stream conditions, all comparisons are made where the upper bound on the body 
angle is restricted so that the sonic line is within the TAT computational 
domain, and the lower bound on the body angle is restricted so that subsonic 
flow exists along the body to the aft corner as required by the present method. 
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Solutions are obtained for the following free-stream conditions: 
Case I: 

Ch = 0.74 
Cfle =0-26 

= 39.09 km/s 
Poo* = X 10-*» kg/m3 
Poo* = 244 N/m2 

Case II: 

Ch = 0.79 

Cfle = 0-21 

VJ» = 28.19 km/s 

P„* = 1.06 X 10-3 kg/m3 

pj* = 0 N/m2 

Case III: 

Ch = 0.79 

CHe = 0.21 

V„» = 31.90 km/s 

Poo* = 7.98 X lO-'* kg/m3 

Poo* = 0 N/m2 

Case IV: 

Ch = 0.79 

CHe = 0.21 

V^» = 39.09 km/s 

P„» = 3.90 X 10-*^ kg/m3 

p » = 0 N/m2 
^ 00 
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Case V: 


C{j = 0.20 

Cfje = 0.80 

V « = 18.13 km/s 
00 


p • = 2.99 X 10“3 kg/m3 
00 

p « = 0 N/m2 
*^00 

(Free-stream pressures equal to zero are acceptable specifications because the 
static pressure contribution to the momentum relationships across the shock are 
extremely small at these velocities.) Execution time for the MIR routine is 
approximately two to three times faster than that for the TAT for all of the 
solutions presented herein. 

Both MIR and TAT use subroutine RAD/EQUIL to calculate the radiative 
heating (ref. 9) . Cases II to V were run with the same version of RAD/EQUIL 
(RAD 69) . An updated version of subroutine RAD/EQUIL (RAD 73) is incorporated 
into the TAT for the purpose of calculating Case I. RAD 73 has been stream- 
lined for faster computations and it treats the hydrogen line contribution to 
the radiation in a more accurate manner. The MIR uses the RAD 69 version for 
all the cases. The input data deck for RAD 69 is modified to reflect some of 
the internal changes in the way RAD 73 treats the hydrogen lines. It is seen 
in figure 7 that both versions of RAD/EQUIL yield radiative heating results 
which agree closely for a given thermodynamic distribution across the shock 
layer typical of Case I. 

Figure 8 presents adiabatic results for surface pressure and shock dis- 
placement distributions corresponding to Case I for a 70° hyperbola with a nose 
radius to base radius ratio of 0.5 (Rji* = 0.223 m) . The MIR prediction of sur- 
face pressure is approximately 7 percent less than predictions made by the more 
exact TAT in the region ahead of the sonic corner. This type of underpredic- 
tion is typical of a single-strip MIR when applied to sonic corner bodies, as 
observed by South (ref. 3) and Jones, et al. (ref. 16) in comparisons made with 
experimental data. Shock standoff distance is also underpredicted in this case, 
especially near the sonic corner. 

The previous case was also computed by including the effect of radiative 
heat transfer. With radiation transfer, the shock layer experiences radiation 
cooling due to the radiation that escapes through the shock. The nonadiabatic 
results cU'e shown in figure 9. Note that there is fair agreement between both 
methods in the stagnation region for shock standoff distance and surface radia- 
tive heating. For this case, the MIR overpredicts the shock standoff distance 
for about one-half of the forebody when compared with the TAT. Also the non- 
adiabatic shock standoff distance is considerably less than the corresponding 
adiabatic value. For example, at the stagnation point, the TAT predicts a 
17 -percent decrease in 6stag while the MIR predicts only an 8-percent 
decrease. At the aft corner of the body, the TAT predicts a 9-percent decrease 
in 6 whereas the MIR predicts only a 3-percent decrease. The comparison of 
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radiative heating predictions closely follows the trend of the shock standoff 
distance predictions. The largest differences between the two prediction tech- 
niques occur along the body near the sonic corner. 

The differences in the percent change of 5 stag brought about by radia- 
tion transfer point out the insensitivity of the MIR relative to that of the 
TAT. This insensitivity is a consequence of the linear approximations in equa- 
tions (4) of properties across the shock layer for a one-strip MIR. Profiles 
of flow field variables as predicted by the two techniques are compared in fig- 
ures 9(b) and 9(c). Note the large gradient in p and h predicted by the 
TAT near the body. On the stagnation line, there is good agreement between p 
and h as predicted by both methods over most of the shock layer. However, 
the radiative cooling near the wall causes the surface values of p and h 
to differ by 19 percent and 42 percent, respectively. ^Percent difference is 

MIR prediction - TAT prediction \ 

defined by x 100.) 

TAT prediction ! 

Pressure predictions agree within 3 percent on the stagnation line. A qua- 
dratic variation of pressure across the shock layer was also tried in several 
cases. Although MIR pressure profiles with a quadratic variation agreed more 
closely with the TAT pressure profiles in the stagnation region, there was 
little effect on the overall flow field predictions. 

Further downstream at a body station s = 1.4, density and enthalpy pro- 
files (fig. 9(c)) follow much the same pattern as in the stagnation region. 

There is a larger percentage difference between the comparisons for density and 
pressure than occurred at the axis of symmetry. The TAT also predicts a larger 
mass flow rate near the body as compared with the MIR solution (fig. 9(c)) 
which is locked into a linear variation across the shock layer according to the 
assumptions made in equations (4). 

When flow field solutions are obtained over less blunt bodies, the MIR 
predictions show larger differences as compared with predictions by TAT. Fig- 
ures 10 to 12 show shock standoff distance and radiative heating distributions 
calculated by both the MIR and the TAT for Cases II, III, and IV over a spheri- 
cally capped cone (0 q = 60°, R(j*/rg = 0.5, and Rij* = 0.223 m) . These condi- 
tions correspond to trajectory points for an entry into a Jovian atmosphere as 
reported in reference 13. 

The comparisons between the two methods are poor. Shock standoff distance 
is accurately predicted in the stagnation region, but it is significantly over- 
predicted by the MIR along the cone in every case. The radiative heating levels 
on the cone are also too high. Some insight into these large discrepancies can 
be gained by looking at the details of the flow profiles for one of these cases. 

Density, enthalpy, and pressure profiles across the shock layer at two 
body locations for Case IV are presented in figure 13. The comparisons here 
between the MIR and TAT are similar to those presented in figure 9. There is 
a large difference between the two methods in the pu profile across the shock 
layer. Because the TAT calculates a larger mass flow rate than the MIR, it 
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follows that the TAT will yield a smaller shock , standoff distance than the MIR 
at this location. In general, the shock layer thicknesses along with the pres- 
sure and temperature levels are critical parameters for determining radiative 
heating levels on planetary probes (ref. 13). The large enthalpy levels near 
the wall and the large shock standoff distances predicted with the MIR combine 
to cause large radiative heating levels compared with TAT results. 

Results of calculations for Case V (spherically capped cone, 6 q = 62°, 
Rff»/rB = 0.5, and Rjj* = 0.223 m) are presented in figure 1^1. This case is one 
for vdiich a two-strip solution could be obtained, as described in the previous 
section. (Free-stream conditions for Case V are those for an entry into Uranus 
as reported in ref. 13.) 

Pressure and temperature distributions are presented in figures 1^(a) 
and l4(b). The two-strip solution improves the MIR pressure prediction as com- 
pared with the TAT, except near the sphere cone junction where a bump occurs in 
the pressure distribution. This bump increased in height as the inner strip 
height decreased. The bump is caused by the discontinuity in body curvature 
at the sphere cone junction. No such irregularities were obtained on any two- 
strip solutions of continuous curvature bodies. Because the irregularity was 
confined to a small area of the flow field, no special differencing or smooth- 
ing procedures were incorporated to try to eliminate this problem. 

The TAT results for temperature variation are shown at nondimensional dis- 
tance of n = 0 (body) and r\ = 0.07. The one-strip MIR agrees closely with 
the TAT prediction at n = 0.07. However, there is a 1000-K temperature dif- 
ference between the TAT results at h = 0 and n = 0.07. This difference is 
indicative of the large gradients in enthalpy profiles near the wall which have 
been observed in all the nonadiabatic results. The two-strip MIR agrees more 
closely with the TAT result for temperature at ri = 0 than with the one-strip 
method. However, even the two-strip method is unable to predict the same large 
decrease in temperature near the wall as the TAT. (See figs. 14(b), 14(e), 
and 14(f).) 

Shock standoff and surface radiative heating distributions are presented 
in figures 14(c) and 14(d), respectively. There is good agreement in the stag- 
nation region among all three methods for shock standoff distance. The two- 
strip MIR and the TAT results agree closely over the entire body. The two- 
strip MIR heating distribution is in good agreement with the TAT results over 
the entire body. The one-strip MIR radiative heating distribution shows fair 
agreement with the TAT in the stagnation region, but it falls below the TAT 
results over the cone. 

The improved comparisons of the two-strip MIR calculations with the TAT 
calculations can be explained by examining the shock layer property profiles 
shown in figures 14(e) and 14(f). There is very little difference between 
the one- and two-strip solutions for the density, enthalpy, pressure, and 
velocity profiles on the stagnation line except near the body. Near the body, 
the enthalpy profile calculated by the two-strip MIR is much closer to the TAT 
result than the profile calculated by the one-strip MIR. However, it is clear 
that even with a two-strip solution, the enthalpy gradient is not nearly as 
great as that calculated by the TAT. Density, enthalpy, pressure, and mass flux 
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profiles calculated by using a two-strip MIR are still in good agreement with 
TAT profiles at a body location s = 1.5 (fig. 1^l(f)). The small bump which 
occurs in the two-strip MIR density profile near the wall is a consequence of 
the linear approximations made on pu and pu^ across the narrow strip. 

^Note that p(q) = £pu(n)]y^[pu2(Ti)^.^ 

Comparisons between the one- and two-strip MIR are presented in figure 15 
for Case V with various cone angles. The relative difference between the one- 
and two-strip radiative heating rate values becomes smaller as the cone angle 
increases . 

The two-strip MIR can resolve the enthalpy and density gradients near the 
wall better than the one-strip MIR. The two-strip MIR also calculates a more 
accurate pressure distribution around the flank of sonic corner bodies. How- 
ever, the two-strip MIR has an extremely limited range of applicability. The 
method can only be used when the sonic line lies near the computational bound- 
ary (fig. 4). In the one comparison case that was available with the TAT, the 
two-strip MIR calculated radiative heating levels that were within 10 percent 
of the TAT values over most of the body. No other cases with entry conditions 
significantly different from the one presented were found that could be calcu- 
lated by using both the TAT and the two-strip MIR. 

One serious drawback of the two-strip method is that the procedure is not 
truly closed. The ability of the two-strip method to converge depends on the 
accuracy of the assumption that the sonic line lies very close to or outside of 
the border of the computational domain at the aft corner of the body. Usually, 
if that assumption is incorrect, the program will fail to converge to any solu- 
tion with the two-strip method. However, some cases have been found where the 
two-strip method converged to a physically unrealistic solution. Also, the two- 
strip method is not always superior to the one-strip method. For example, there 
is little difference in the radiative heating results between the one- and two- 
strip methods for large cone angles (fig. 15). Consequently, the range of con- 
ditions over which the two-strip method can offer any significant improvement 
over the one-strip method is small. 


CONCLUDING REMARKS 

The one-strip method of integral relations is unreliable in making radia- 
tive heating calculations over planetary probes with sonic corners. The method 
yields surface pressures which are too low near the sonic corner and is unable 
to resolve enthalpy and density gradients near the wall. As the solution is 
marched downstream from the stagnation line, these errors combine to produce 
gross distortions of the radiative heating levels to the surface. 

A modification of a two-strip method of integral relations, as described 
herein, does a better job of resolving the enthalpy and density gradients near 
the wall than the one-strip method; however, the gradients predicted by the 
two-strip method are still not quite as severe as those predicted by the time 
asymptotic technique. The two-strip method does not distort the pressure dis- 
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tribution around the body as severely as the one-strip method. In the one 
detailed comparison case with the time asymptotic technique, the two-strip 
method of integral relations accurately predicted the radiative heating levels 
over the entire body. 

Very little difference was observed in the radiative heating distributions 
between the one- and two-strip methods for very large cone angles (0 q > 75°) . 
Thus, for the very blunt bodies, there seems to be little advantage in going to 
the two-strip method. There is also the possibility that the two-strip method 
can converge to a physically unrealistic solution because the procedure is not 
truly closed. 

Consequently, the one-strip method of integral relations cannot be relied 
on to obtain accurate estimates of radiative heating over planetary probe 
vehicles for the conditions described herein. In some instances, a two-strip 
method can be used to improve these estimates. However, the small range over 
which the two-strip method can offer significant improvement over the one-strip 
method makes it unrealistic to pursue further development of this technique. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 21, 1977 
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APPENDIX 


LOCALLY VALID RELATIONS FOR EQUATION OF STATE AND SPEED OF SOUND 

Subroutine RAD/EQUIL (ref. 9) is used to obtain density, the mole fraction 
of the considered species, the ratio of specific heats, and temperature as a 
function of pressure and static enthalpy for arbitrary gas mixtures. An equa- 
tion of state for a real gas in chemical equilibrium must be coupled with the 
other governing equations to obtain a flow field solution. A successive sub- 
stitution algorithm is time consuming, and in the region near the sonic line 
singularity, the successive substitution can diverge. It is desirable to 
derive a local equation of state at each point in the flow field which is 
valid over the small range of suiticipated pressures and enthalpies at that 
point. This analytic equation of state is then coupled with the other govern- 
ing equations and the entire system of equations is solved by using the qua- 
dratically convergent Newton-Raphson technique. It is assumed that a locally 
valid equation of state can be written as 

p = Ap^hC (A1) 

Choose Pg, pjj, and Pq equal to three values of pressure within the 
anticipated pressure range. Choose h^, h^, and hg equal to three values 

of enthalpy within the anticipated enthalpy range. (The anticipated pressure 
range near the axis of symmetry is bounded by the stagnation pressure and the 
pressure behind the shock vrtiich can be obtained directly from RAD/EQUIL. The 
anticipated enthalpy range near the axis of symmetry is bounded by the stagna- 
tion enthalpy and static enthalpy behind the shock. Anticipated ranges of 
pressure and enthalpy at subsequent computational points are extrapolated from 
converged values of pressure and enthalpy upstream of the points.) 

Let 

Pa = P(Pa»ha) 

Pb = PCPb'^’b) 


Pc ” P(Pc»^C^ 


where p(p,h) is obtained from subroutine RAD/EQUIL. If equation (A1) is to 
be valid, the following equations must be true: 


Pa = APa®ha^ 

Pb = APb^hb^ > 
Pc ~ ApQ®hQ^ 


(A2) 
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Taking the logarithm of both sides of equations (A2) euid solving for A, 
B, and C yields 


[loge (Pa/Pc) loge (Pa^Pb) * l°ge (Pa/Pb^ l°ge (Pa/Pc)] 

C = ; (A3) 

[loge (pa/Pc) loge ^^a/hb) - logg (Pa/Pb) loge (ha/hc)] 

[loge (ha/hc) logg (Pa/Pb) - loge (ha/hb) logg (Pa/Pc)] 

B = (a4) 

[loge (ha/hc) loge (Pa/Pb> - loge (^a/hb) loge (Pa/Pc)] 

A = PaPa^^^^a'^ (A5) 

Values for A, B, and C are calculated in this manner at the boundary 
of every strip and at every integration step around the body. These values are 
updated at every pass around the body. (Because the pressure and enthalpy lev- 
els in the stagnation region do not substantially change during the last few 
iterations, it should be possible to save values of A, B, and C, and thus 
save additional calls to RAD/EQUIL. However, the present program was never 
optimized in this manner.) 

A comparison of results for density obtained by using a local equation of 
state and a direct use of RAD/EQUIL are presented in table I. 

It is necessary in this analysis to be able to locate the sonic line in 
the flow field. A speed of sound must therefore be calculated at every compu- 
tational point in the flow field. The speed of sound is defined as 

a» = 

However, y> Rg*> and T* are functions of pressure and enthalpy for arbi- 
trary gas mixtures . Since a local equation of state must be determined at 
every computational point , it is convenient to determine a local equation for 
the speed of sound simultaneously, and thus eliminate an additional call to 
RAD/EQUIL . 

It is again assumed that a local equation for the speed of sound can be 
written as 

a = DpRpF 

By following the development in the previous section, it is found that 


[loge (Pa/Pc) loge (aa/^b) " logg (Pa/Pb^ loge (aa/ac^l 

F = : = (A6) 

[loge (Pa/Pc) loge (Pa/Pb) ~ loge (Pa/Pb^ loge (Pa^Pc)] 
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floge (Pa/Pc) ^°8e (^a/ab) - logg (Pa^Pb) l°ge (ag/ac)] 

E (A7) 

[logg (Pa/Pc) loge (Pa/Pb) - logg (Pa/Pb^ loge (Pa^Pc)] 


D = aaPa"%>a‘^ (A®) 

where 

aa “ aCpg^jhg) 
ab “ a(p^j,h||j) 

Bq - s(Pq jhg) 

and p's, p's, and h's with subscripts a, b, and c are defined in the 
previous section. A comparison of results for speed of sound -obtained by use 
of a local equation for the speed of sound and a direct use of RAD/EQUIL are 
presented in table I. 
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TABLE I.- COMPARISON OF APPROXIMATE RESULTS CALCULATED BY LOCALLY VALID 


RELATIONS WITH EXACT RESULTS CALCULATED BY RAD/EQUIL 


3 

n 

Density 

Speed of sound 



Approximate 

Exact 

Percent 

Approximate 

Exact 

Percent 





difference 



difference 

Adiabatic 

0 

0 

10.863889 

10.863889 

0.00000000 

0.32746270 

0.32746270 

0.00000000 


1 

10.413217 

■ 

10.415360 

.02057538 

.32580473 

.32577076 

.01042760 

.5 

0 

10.498678 

10.498678 

.00000000 

.32619309 

.32619307 

.00000610 


1 

10.371731 

10.371735 

.00003857 

.32341920 

.32341908 

.00003710 

1.0 

0 

10.151106 

10.151106 

. 00000000 

.32488561 

.32488561 

.00000000 


1 

10.294323 

10.294323 

.00000000 

.31950323 

.31950322 

.00000310 

1.5 

0 

9.669698 

9.669696 

.00002070 

.32299906 

.32299910 

.00001238 


1 

10.185146 

10.185141 

.00004910 

.31421897 

.31421911 

.00004455 

2.0 

0 

8.200939 

8.200835 

.00126820 

.31672379 

.31672656 

.00087457 


1 

9.941612 

9.941542 

.00070410 

.30319808 

.30320109 

.00099274 

Nonadiabatic 

0 

0 

12.977966 

12.977966 

0.00000000 

0.29983841 

0.29983841 

0.00000000 


1 

10.413217 

10.415360 

.02057538 

.32580473 

.32577077 

.01042450 

.5 

0 

13.212973 
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Figure 2.- Flow chart of radiative heating coupling procedure. 
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Figure 5.- Comparison of one- and two-strip MIR results for shock shape over 
70° spherically capped cone with = 0*5, M = 6.02, and y - 1.67. 
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Figure 6.- Initial conditions yielding sonic velocity on corner of body for 
two-strip solution with M = 6.02, y = 1.67, and 9q = 7 



RAD 69 (modified) 



Figure 7.- Comparison of shock layer results using the RAD 73 and the modified 
RAD 69 radiation subroutines for identical pressure and enthalpy profile 




(a) Surface pressure distributions. 


Figure 8 



(b) Shock standoff distributions. 

- Comparison of MIR calculations with TAT results. Case I; adiabatic 











(a) Shock standoff distributions. 



(b) Radiative heat-transfer distributions. 

Figure 10.- Comparison of nonadiabatic results for Case II 

Rjj/rB = 0.5; and Rjj* = 0.223 n. 



(a) Shock standoff distributions. 
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(b) Radiative heat-transfer distributions. 


Figure 11.- Comparison of nonadiabatic results for Case III 

and Rm* = 0.223 m. 



(a) Shock standoff distributions. 
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(b) Radiative heat-transfer distributions. 

Figure 12.- Comparison of nonadiabatic results for Case IV 

Rjg/rB = 0.5; and Rfj* = 0.223 
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(a) Pressure distribution. 



(b)' Temperature distribution. 

Figure 14.- Comparison of nonadiabatic results for Case V. 

%/i^B = •^•5; and Rm» = 0.223 m. 
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(e) Shock layer profile properties at s = 0. 
Figure 14.- Continued. 
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Figure 15.- Effect of cone angle on radiative heat transfer for case V 

Rjg/rB = 0.5; and Rjg* = 0.223 ni. 
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